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Abstract. We report on the infrared limit of the quenched lattice Landau gauge gluon and ghost propaga- 
tors as well as the strong coupling constant computed from large asymmetric lattices. The infrared lattice 
propagators are compared with the pure power law solutions from Dyson- Schwinger equations (DSE). For 
the gluon propagator, the lattice data is compatible with the DSE solution. The preferred measured gluon 
exponent being ~ 0.52, favouring a vanishing propagator at zero momentum. The lattice ghost propagator 
shows finite volume effects and, for the volumes considered, the propagator does not follow a pure power 
law. Furthermore, the strong coupling constant is computed and its infrared behaviour investigated. 

PACS. 12.38.-t Quantum chromodynamics - 11. 15. Ha Lattice gauge theory - 12.38.Gc Lattice QCD 
calculations - 12. 38. Aw General properties of QCD - 14.70.Dj Gluons - 14.80.-j Other particles 



1 Introduction 

In the pure gauge SU(3) Yang-Mills theory, a number 
of authors has been using first principles approaches, i.e. 
Dyson-Schwinger equations (DSE) and lattice QCD meth- 
ods, to investigate the infrared gluon and ghost propaga- 
tors in Landau gauge, respectively. 
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and the strong coupling constant [T] defined as 

as{k^) = as{fi^)Z'^f^^^t{k^)Zgiuonik'^); (3) 

Zghostik^) = fc2G(fc2) and Zgi„o„(fc2) = k^D{k^) are the 
ghost and gluon dressing functions. In part, these stud- 
ies have been trigerred by the solution of the DSE [2] 
which assuming infrared ghost dominance and infrared 
finiteness of the loop-integrals predicts pure power laws 
for the dressing functions 
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Moreover, the DSE solution relates the exponents of the 
two propagators, k' = 2k, and predicts a finite strong 
coupling constant at zero momentum. The infrared solu- 
tion predicts a vanishing gluon propagator and an infi- 
nite ghost propagator at zero momentum (k > 0.5). The 
infrared behaviour of the propagators can be related to 
gluon confinement criteria [3^. Looking at this particular 
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Fig. 1. DSE gluon and ghost dressing functions versus the 
pure power law solutions. The data is taken from [4]. 



DSE solution for the gluon propagator, the comparison 
with the pure power law, see figure [H shows that the pure 
power law is valid only for momenta k < 200 MeV. Notice 
that for the ghost, the deviations from a pure power law 
start earlier. 

One should have in mind that the solution discussed 
above is a particular solution of the DSE. Indeed, there 
are, in the literature, different types of solutions for the 
DSE fF. Moreover, the authors of 6J have investigated a 
generalization of the conditions assumed in [2j and found 
alternative behaviours for the infrared gluon and ghost 
propagators. Given the different scenarios, it would be im- 
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portant if lattice QCD can provide additional input, help- 
ing to check if any of the proposed solutions reproduces 
the lattice data. For a recent discussion on the infrared 
behaviour of the gluon and ghost propagators see [7] . 

For the lattice it is a challenge to investigate such 
low momenta. A possible way out is to consider large 
asymmetric lattices, which allows to access the momenta 
required in such an investigation [8l[9lfT0l[TT] . Of course, 
there are lattice effects that have to be carefully estimated. 
Here we report on the gluon, ghost and strong coupling 
constant computed from large asymmetric lattices. 



2 Gluon Propagator 

In [11] we have investigated the infrared lattice Landau 
gauge gluon propagator. For notation and definitions see 
the above cited article. The simulation uses the Wilson 
action with (3 = 6.0 and combines a number of asymmet- 
ric lattices x 256, L = 8, 10, 12, 14, 16, 18 to allow for 
L +O0 extrapolation. In what concerns the time direc- 
tion, the 16*^ X 256 and 16"^ x 128 gluon propagators were 
compared and no deviations were observed. This suggests 
that we have a sufficient number of points in the time 
direction. 

For the lattices with the largest time direction, the 
minimum momentum being accessed is 47 MeV, while for 
the 16"^ X 128 the minimum momentum is 98 MeV. For 
the largest lattices with T = 256 a pure power law in the 
infrared is observed. In contrast, for the lattice of 16^^ x 
128, power law fits are very poor with ^/d.o.f > 10. 
This result is not surprising, given the smallest range of 
momenta available for this lattice (97 - 294 MeV) and the 
validity of the pure power law. Nevertheless, the 16"^ x 
128 will allow us to estimate Gribov copies effects, with 
reasonable computational effort. 

Let us summarize the results reported in [llj using 
the simulations with the largest time extension. Asym- 
metric lattices display finite size effects (see figures 4 and 
5 of [H]). However, the extrapolation towards L +oo 
is smooth. The extrapolated gluon propagator was well 
reproduced by a pure power law, with measured k = 
0.49 — 0.52. Clearly, the lattice data favors k ~ 0.52, 
in agreement with other theoretical estimates. Unfortu- 
nately, one cannot yet provide a definitive answer concern- 
ing the value of the gluon propagator at zero momentum. 
We are currently engaged in trying to give such an answer. 

For the 16'^ x 128 lattice, Gribov copies effects were 
estimated by comparing the gluon propagator computed 
from 164 configurations gauge fixed as in [11] (SD in the 
figures) and the same gauge configurations gauge fixed 
using the method described in pjj (CEASD in the fig- 
ures), which aims to estimate the absolute maximum of 
the gauge fixing function. The gluon propagator given by 
the two methods is, within the statistical precision of the 
simulation, the same (see figure [2]), with the various fits 
to the lattice data (infrared, ultraviolet, all momenta) re- 
producing similar results. Therefore, no visible effects of 
Gribov copies are observed on the gluon propagator. 
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Fig. 2. Bare gluon propagator for 16^ x 128. 

3 Ghost Propagator 

The ghost propagator [13] was computed using two differ- 
ent methods [14,..15j. The first allows, by solving a linear 
system, to access all momenta. Its drawback being that, 
by computing G{x,0), the statistical errors in the prop- 
agator are much larger. The second method requires an 
inversion of the same linear system per momentum, which 
is more computationally demanding. However, for the mo- 
menta considered, it determines the propagator with bet- 
ter statistical accuracy. Computing the propagator using 
two methods provides a valuable cross-check. 

In the following, for the ghost propagator the colour 
average was always performed. For the propagator com- 
puted with the [13] method, in order to reduce the sta- 
tistical error, seven different sources for the linear system 
were considered and their result averaged. 

For the ghost propagator, one observes sizeable finite 
size effects (see figure in agreement with the SU(2) 
study [IB] , and a very sign of clear Gribov copies effects 
for a large range of momenta. 

In what concerns finite volume effects, figure |3] shows 
that the ghost dressing function is enhanced when the 
lattice volume is increased. Moreover, the larger volumes 
have a larger dressing function and, for the smallest mo- 
menta, the derivatives of the dressing function seem to 
become smaller when the volume increases. The effect of 
Gribov copies, seen in figure jS] cannot be eliminated by 
the renormalization of the propagator. 

To study the compatibility of the ghost propagator 
with the DSE solutions, the lattice data was fitted to pure 
power laws and naive corrections to the pure power law. 
It turns out that the data are not described by any of the 
functions considered. Given the results of a similar study 
for the 16"^ x 128 gluon propagator, this result is not a 
complete surprise. 



4 Strong Coupling Constant 

The strong coupling constant as defined by equation ([3]) is 
displayed in figure 11 for the 16^ x 128 lattice and CEASD 
gauge fixing method. For the other simulations and the 
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Fig. 3. Bare gliost dressing functions for a 16'^ x 128 and 
X 256, L = 10, 12 lattices. 
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Fig. 5. Small momenta strong coupling for all the simulations. 
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Fig. 4. Strong couphng constant for 16^ x 128 configurations 
for the CEASD gauge fixing method; "p2p" ( "pws" ) stands for 
the ghost computed using [M] ([15j) method. 



SD method, the measured coupling constants are, quali- 
tatively, similar. Figure [5] show the strong coupling con- 
stant for small momenta, for all the simulations considered 
previously. 

The simulations show a decreasing lattice strong cou- 
pling constant for small momenta, in agreement with pre- 
vious lattice calculations 17,18 . From our simulation, it 
is not clear if as(0) is vanishing. We have tried a number 
of fits to the strong coupling constant as a function of the 
momentum, some requiring as(0) = and some leaving 
as(0) as a free parameter, and the data can be equally de- 
scribed by the two situations. For small momenta, only for 
the 16^ X 128 data and the CEASD gauge fixing algorithm, 
is as described by a pure power law as(q^) = {q'^)^" , with 
Ka = 0.69 suggesting a vanishing strong coupling constant 
at zero momentum. However, notice that the data show 
that as increases with lattice volume as seen in figure 
[5] and is affected by Gribov copies. Further studies will 
hopefully lead to definitive conclusions. 
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